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Abstract 

Charged hadron differential distributions from muon-induced deep 
inelastic scattering, DIS, on a 6 LiD target are presented as function 
of the DIS variables XBj, Q 2 , W 2 and the hadron variables px and z. 
They can be used as benchmark to verify the reliability of theoretical 
model. The p T distributions are fitted with a Gaussian function at 
different kinematic intervals. With a Gaussian ansatz for the trans- 
verse momentum dependent parton distributions, TMDs, the intrinsic 
transverse momentum of the partons is extracted. 

1 Introduction 

Semi Inclusive Deep Inelastic Scattering (SIDIS) reactions provide 
much information about the structure of the nucleon and the hadroniza- 
tion of partons. Parton intrinsic transverse momentum affects the 
transverse momentum of the produced hadron. The kinematic vari- 
ables of the produced hadron are shown in figured! where the trans- 
verse momentum is determined w.r.t. the virtual photon. The trans- 
verse momentum distributions from unpolarized scattering are often 
the first verification of the reliability of a theoretical model. The shape 
of the p T distributions depends on many effects, e.g. contribution from 
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Figure 1: Kinematic variables of virtual photon interaction with a parton 
with intrinsic transverse momentum k± and its hadronization. The transverse 
momentum of the observed hadron, p?, is defined with respect to the virtual 
photon. 

intrinsic transverse momentum k± to the TMDs and transverse mo- 
mentum of hadron p± to the fragmentation functions (FFs) and gluon 
radiation. These effects depend on kinematic variables such as the 
Bjorken variable XBj, the invariant mass squared W 2 of the hadronic 
system, the negative 4-momentum squared Q 2 of the virtual photon 
and the virtual photon energy fraction z carried by the hadron. Af- 
ter integrating over the azimuthal angles, four variables are needed to 
describe the kinematics of the measured hadron: two inclusive (Q 2 , 
xsj) and two hadronic (pt, z). 

2 The COMPASS Experiment 

The COMPASS experiment has been set up at the M2 muon beam 
line of the CERN SPS [I]. Polarized 160 GeV muons with an inten- 
sity of 2 • 10 8 /x/4.8s spill and a polarization of 80% are scattered off 
a polarized 6 LiD target. The target consists of two cells of opposite 
polarization which was reversed every 8 h. The unpolarized sam- 
ple is therefore the combination of the data from the two cells. The 
COMPASS spectrometer is a large acceptance two-stage spectrometer 
which covers the kinematic range from quasi-real photo production 
to the DIS region. Both stages use hadron calorimeters and absorber 
walls for muon identification. The data presented here were taken 
during the year of 2004. 
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3 Hadron kinematic distributions 

The charged hadron identification is kept as as simple as possible. 
The particles coming out of the primary vertex are either identified 
as hadron or muons. Prom these hadrons, the selection requires that 
they create signals in the detector situated upstream and downstream 
of the first magnet. This ensures that the track momentum and charge 
are well defined by the bending of the magnetic field. The COMPASS 
ability to identify hadrons, with a RICH detector, was not used, but 
is intended for further analysis. 

In order to correct for event losses caused by the non-uniform accep- 
tance of the COMPASS spectrometer, a full Monte Carlo (MC) simu- 
lation has been performed. The events were generated with LEPTO, 
transported through the COMPASS detector simulation program COMGEANT 
and the reconstruction software CORAL. From this MC sample, 4- 
dimensional acceptance tables have been determined. Although very 
similar, positive and negative hadrons have different tables. The sys- 
tematic error has been estimated to 5%. Only statistical errors are 
shown in the figures. In this analysis, the hadrons are separated into 
23 intervals in Q 2 (from 1 to 10 GeV/c 2 ) and x Bj (from 0.004 to 0.12) 
further subdivided into 8 intervals in z (from 0.2 to 0.8). 

4 Results 

Hadron muoproduction has been studied for many years and the EMC 
experiment [5] covered a similar kinematic range as COMPASS. An 
interesting comparison with previous data is the ratio of positive and 
negative hadrons because the acceptance is canceled to a good ap- 
proximation. The hadron multiplicity ratios are shown in figure [2j 
COMPASS results show clearly the z and £ ^-dependence, where the 
fraction of positive hadrons increases with XBj (getting closer to the 
valence region) and z (more related to the struck parton). This agrees 
with the model of valence quarks where the positive quarks have a 
higher electric charge. 

4.1 "p\ distributions 

According to [I], the average over all -p\ should depend linearly on the 
center of mass energy squared, s. They have verified their prediction 
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Figure 2: Charged hadron multiplicity ratios dn + /dn as function of z for 
EMC [5J and COMPASS. 

with results from three fixed target experiments: JLab, HERMES 
and COMPASS. The value used for COMPASS was not corrected for 
acceptance, the value is corrected here as shown in figure El The 
authors of [I] note that the average °p\ should depend linearly on W 2 
rather than s. The dependence is shown in figure [JJ which is more 
compatible with a linear dependence on W than on W 2 . The relation 
is not well established and, as mentioned in [JJ, the linear dependence 
on s for Drell-Yan which inspired their SIDIS prediction, could also 
be a linear dependence on ^/s. Figure H] suggests the latter is more 
accurate. 

4.2 Gaussian fit of the "p\ distributions and in- 
trinsic transverse momentum 

Using the Gaussian ansatz which assumes a Gaussian distribution of 
the intrinsic transverse momentum and of the transverse momentum 
acquired during fragmentation, the cross section is proportional to (cf. 

my- 

, 2 , ,2 OcX>g/g(sBj)A,(z) , (1) 



dx Bj dQ 2 dzdp 2 T ^ qJqy Jl qy ' tt{ p 2 
where 

(p 2 T ) = (p 2 ± ) + z 2 (k 2 ± ). (2) 

The functions f q (xBj) and D^(z) are the usual (integrated) distribu- 
tion and fragmentation functions, respectively. The distributions and 
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Figure 3: The left figure shows the differential p\ distribution of hadrons 
with 0.5 < z < 0.6. It is used to determine the corrected COMPASS averaged 
over all p\. This value is to correct the figure on the right from [3], where 
three experiments show the s-dependence of the average p\. The red star 
COMPASS point is the value from this analysis, the black dot, COMPASS 
(Schw.), is the value inferred by jl] from data uncorrected for acceptance. 
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Figure 4: Charged hadron average over all p\ vs W 2 . Note that contrary to 
all the following figures of this article, (pj) defines a standard average over 
all p\. 
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Figure 5: "p\ distributions fitted by a Gaussian for (1 < Q 2 < 1.5 (GeV/c) 2 , 
0.006 < x Bj < 0.008) and (6 < Q 2 < 10 (GeV/c) 2 , 0.07 < x Bj < 0.12) 
subdivided into 8 z intervals. 

fits for two (Q 2 , XBj) intervals and for all z intervals are shown in 
figure[5j The fit is performed on the low pt interval [0.1, 0.85] GeV/c 
in order to stay away from pQCD effect. The fitted (p 2 -) as function 
of XBj are shown in figure [6] for low and high z. A similar behavior 
was already observed by HERMES in [6j for the average pr (not from 
a fit but from a standard average). It is interesting to compare the 
average p\ from the previous section (figure H]) with the fitted {p 2 -) for 
the middle z shown in figure [71 Contrary to the average aver all pj,, 
there is no clear VF 2 -dependence of the fitted {p 2 -) (for 0.1< px <0.85 
GeV/c) which is suppose to be unaffected by pQCD. 

The z 2 -dependence of the fitted (pj) is of particular interest be- 
cause of its relation to the intrinsic transverse momentum through 
equation fl2J). The fitted {pj) for different z 2 for two (Q 2 , xsj) in- 
tervals are shown in figure [8j The relation between (pj) and z 2 is 
certainly not linear as in equation ([2]). If a z-dependence of the trans- 
verse momentum acquired during fragmentation, p±, is added such 
that 

(p 2 T )=z a (l-zf(p 2 ± )+z 2 (k 2 ± ), (3) 

where a = 0.5 and /3 = 1.5, the relation can be nicely fitted as shown 
in figure [HJ From this fit, the intrinsic momentum, (k 2 ^), for various 
(Q 2 , XBj) intervals can be extracted. The extracted (kj_) as function 
of Q 2 are shown in figure [9] and as function of XBj, comparing results 
from positive and negative hadrons, are shown in 1101 
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Figure 6: Fitted (p?) vs xbj for different Q 2 intervais for low (left column) 
and high (right column) z. In the present figure, the (p 2 -) results from a fit 
over p T < 0.85 GeV/c. 
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Figure 7: Fitted vs W 2 at medium z ([0.5, 0.6]) for different Q 2 intervals. 
This is to be compared with figure H] where the average over all p\ is plotted. 
In the present figure, the (p^) results from a fit over px < 0.85 GeV/c. 
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Figure 8: Fitted (p?) vs z 2 for two (Q 2 , %Bj) intervals. The fit function is 
given by equation ([3]). The dotted green line is the result of a fit, performed 
by [3], to data from many experiments. In the present figure, the (p 2 ^) results 
from a fit over px < 0.85 GeV/c. 
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Figure 9: Extracted (k\) vs Q 2 for various xsj intervals. The dotted green 
line is the result of a fit, performed by j3], to data from many experiments. 
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Figure 10: Extracted (k±) vs xbj for positive and negative hadron compar- 
ison. The dotted green line is the result of a fit, performed by [3], to data 
from many experiments. 
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5 Conclusion 



The differential p\ distributions of charged hadrons produced by muons 
scattered off a 6 LiD target have been determined for various kinematic 
intervals. The low p\ have been fitted with a Gaussian at different z 
such that the intrinsic transverse momentum could be extracted in the 
framework of the Gaussian ansatz. The non linear relation between 
the fitted {p\) and z 2 have been reproduced by adding a z-dependence 
of the transverse momentum acquired during fragmentation. The ex- 
tracted (k±) shows a clear dependence on Q 2 and a possible depen- 
dence on XBj, although less conspicuous for higher Q 2 , where the 
framework is more reliable. Also, is systematically higher for 

positive hadrons compared to negative hadrons. This suggests a fla- 
vor dependence of the intrinsic transverse momentum. This behav- 
ior could be further investigated using COMPASS ability to identify 
hadrons; kaon identification could provide access to characteristics of 
the strange quark TMDs. 
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